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Abstract. Connections between microscopic follow-the-leader and macroscopic 
fluid-dynamics traffic flow models are already well understood in the case of 
vehicles moving on a single road. Analogous connections in the case of road 
networks are instead lacking. This is probably due to the fact that macroscopic 
traffic models on networks are in general ill-posed, since the conservation of 
the mass is not sufficient alone to characterize a unique solution at junctions. 
This ambiguity makes more difficult to find the right limit of the microscopic 
model, which, in turn, can be defined in different ways near the junctions. In 
this paper we show that a natural extension of the first-order follow-the-leader 
model on networks corresponds, as the number of vehicles tends to infinity, to 
the LWR-based multi-path model introduced in [4, 5]. 


1. Introduction. Traffic flow can be described at different scales, depending on 
the level of details one wants to observe. Typically, three scales of observation 
can be adopted: microscopic (single vehicles are tracked), mesoscopic (averaged 
quantities like density and mean velocity are tracked but car-to-car interactions are 
not completely lost) and macroscopic (only averaged quantities are observed). 

Limiting our attention to differential models, the most famous microscopic mod¬ 
els are those of follow-the-leader type, also known as car-following models [20, 27]. 
In such a models, the dynamics of each vehicle depend on the vehicle in front of it, so 
that, in a cascade, the whole traffic flow is determined by the dynamics of the very 
first vehicle (the leader). The oldest macroscopic model is instead the LWR model 
[25, 28], which dates back to 1955 and is inspired by the laws of fluid-dynamics. 

Generally speaking, microscopic models are considered more justifiable because 
the behavior of every single vehicle can be described with high precision and it 
is immediately clear which kind of interactions are considered. On the contrary, 
macroscopic models are based on assumptions which are hardly correct or at least 
verifiable (among others, the actual validity of the continuum hypothesis and the 
negligibleness of car-to-car interactions). As a consequence, it is often desirable 
establishing a connection between microscopic and macroscopic models in order to 
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justify and validate the latter on the basis of the verifiable modeling assumptions 
of the former. 

Connections between microscopic follow-the-leader and macroscopic fluid-dynamics 
traffic flow models are already well understood in the case of vehicles moving on a 
single road. Analogous connections in the case of road networks are instead lacking. 
This is probably due to the fact that macroscopic traffic models on networks are 
in general ill-posed, since the conservation of the mass is not sufficient alone to 
characterize a unique solution at junctions. This ambiguity makes more difficult 
to find the right limit of the microscopic model, which, in turn, can be defined in 
different ways near the junctions. 

Goal. In this paper we propose a very natural extension of a first-order follow- 
the-leader model on road networks and then we prove that its solution tends to the 
solution of the LWR-based multi-path model introduced in [4, 5] in the limit, i.e. 
as the number of vehicles tends to inhnity while their total length is kept constant. 
The limit is proved extending to networks the results already existing for a single 
road, and it is then confirmed by numerical experiments. It is useful to recall 
that the multi-path method is able to select automatically an admissible solution 
at junction, thus resolving ill-posedness issues. We also recall that the solution 
selected by the multi-path method does not match the one obtained by maximizing 
the flux at junction. Therefore, the connection with the microscopic model promotes 
the solution computed by the multi-path method as “more natural”, while the one 
computed by maximizing the flux should be seen as the “most desirable”, to be 
achieved by means of ad hoe traffic regulations. 

Relevant literature. The literature about microscopic and macroscopic traffic 
flow models is huge and a detailed review is out of the scope of the paper. For 
a quick introduction to the held we suggest the book [19] and the surveys [2, 20]. 
Regarding hrst-order models on a single unidirectional road, the micro-to-macro 
limit was already deeply investigated by means of different techniques: the papers 
[7, 29] use standard techniques coming from the theory of conservation laws; the 
paper [9] instead proves the limit relying also on measure theory. The microscopic 
solution is interpreted as an empirical measure which is proven to converge to the 
entropy solution of the macroscopic model in the I-Wesserstein topology; hnally, 
the papers [8, 13] attack the problem exploiting the link between conservation laws 
and Hamilton-Jacobi equations. 

Micro-to-macro limit for second-order models was instead investigated in [1, 18], 
where the Aw-Rascle model is derived as the limit of a second-order follow-the- 
leader model. 

Macroscopic-only traffic models on networks were deeply investigated starting 
from [24]. A complete introduction can be found in the book [16], which discusses 
several methods to characterize a unique solution at junctions. Let us also men¬ 
tion the source-destination model introduced in [15] (see also [21]) and the buffer 
models [14, 17, 23]. Recently, the LWR-based multi-path model on networks was 
introduced in the paper [4], together with a Godunov-based numerical scheme to 
solve the associated system of conservation laws with discontinuous flux. The re¬ 
lationship between the multi-path model and more standard methods (like, e.g., 
the one proposed in [6] based on the maximization of the flux at junction) was 
investigated in [5]. 

To our knowledge, there is no systematic theory about the extension of the follow- 
the-leader models on networks. It is plain that at the microscopic level one can easily 
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reach a high level of detail, including junctions with spatial extension (non point), 
multi-lane roads, multi-class vehicles, traffic lights and priorities. Several highly 
sophisticated simulators are available since many years (free and commercial), see, 
e.g., [10, 30] and references therein to have an idea of the models and methods 
commonly used. Nevertheless, it is unclear which average flux is actually observed 
at junction by the many-particle limit of any car-following model. 

Beside this, let us also mention that the relationship between microscopic and 
macroscopic models was exploited to create hybrid models, see, e.g., [26]. In such 
a models the averaged quantities are observed where a detailed description is not 
needed (e.g., far from the junctions) while microscopic dynamics are considered 
elsewhere. However, this approach gives no clue about the macroscopic behavior of 
the microscopic model at junctions. 

Paper organization. In section 2 we introduce the follow-the-leader model and 
the LWR model on a single road. We also recall the existing results about the 
micro-to-macro limit on a single road following [7, 9, 29]. In section 3 we extend the 
follow-the-leader model to networks, and in section 4, which is the core of the paper, 
we show the relationship between the follow-the-leader model on networks and the 
LWR-based multi-path model. In section 5 we present fully discrete algorithms for 
the numerical solution to the equations related to the models previously discussed 
and finally in section 6 we confirm our findings by means of some numerical tests. 

2. Background and previous results. Let us describe the first-order follow-the- 
leader model we will use as main ingredient in the rest of the paper. We assume 
here that vehicles move on a single infinite road with a single lane. Vehicles are 
initially located one after the other and cannot overtake each other. We denote by 
n € N the number of vehicles, by £„ > 0 the length of the single vehicles, and by £ 
their total length. We have 



This relation is crucial for the micro-to-macro limit since it translates the fact that 
the total length does not change when the number of vehicles tends to infinity, 
because the cars “shrink” accordingly. 

We denote by yi(t) the position of the i-th car and we assume that at the initial 
time t = 0 cars are labeled in order, i.e. 2/i(0) < 2 / 2 ( 0 ) < ... < 2/n(0). This 
guarantees that the {i + l)-th car is just in front of the i-th one. Moreover, it is 
assumed that at t = 0 cars do not overlap, i.e. 2/i+i(0) — 2/i(0) > in, i = 1, ■ ■ ■ ,n — 1. 
We are now ready to introduce the model, described by the following system of 


DDEs 


y,{t) = w{Si{t)), i = l,...,n-l 
iiniji) = Umax; 


( 2 ) 


where 


5i{t) := yi+i{t) - yi{t) 


and w is such that 



maxj 7 



(3) 


with V G (^^([0,1]; [0, Umax]) any function such that 

u'(r) < 0, u(0)=Umax, u(l) = 0. 


'max. 


( 4 ) 
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The n-th car is the leader and it is assumed to move at maximal velocity Vmux > 0. 
Note that the properties introduced above guarantee that the cars do not overlap 
at any later time t > 0, see [9, Lemma 1]. 

The macroscopic limit of the previous model is given by the well known LWR 
model, which describes the evolution of the average (normalized) density of vehicles 
p(t, x) : ([0, +CX)) X R) —>• [0,1] by means of the following conservation law 


dtp + dx{pv{p)) = Q, t > 0, X G R. (5) 

Relationship (3) makes the link between the two models. 

In order to recall precisely the results about the correspondence between the two 
models, we need to introduce first the natural spaces for the macroscopic density p 
and for the vectors of vehicles’ positions y = (j/i,..., ?/„), at any fixed time: 

R := |r G L^(R;[0,1]) : L r{x)dx = C and supp(r) is compact 


and 


Yn:={y€ 




y^+l Ui ^ Vz — 1 , . . . , 72 1 




We also introduce the operators En 
respectively as 

yn = max(supp(r)). 


R ^ Yn. and 




R, defined 


En 


and 


•)] :=y = 


Pi = max < z € 


rVi+i 


r(x)dx = in 


n-l 


Cn[y] fn — ^ . X[yi,yi+i)^ 

i=l 


i = n — 1,..., 2,1 

( 6 ) 

( 7 ) 


where x —> x/(x) is the indicator function of any subset / C R. The discretiza¬ 
tion operator En acts on a macroscopic density p(t, •), providing a vector of posi¬ 
tions y(t) = En[p{t, •)] whose components partition the support of the density into 
segments on which p has fixed integral ^n■ On the contrary, the operator Cn an¬ 
tidiscretizes a microscopic vector of positions y, giving a piecewise constant density 

Pn{t,-) = Cn[y{t)]. 

We are now ready to state the main result about the convergence of the micro¬ 
scopic model (2) to the macroscopic one (5). The proof can be found in [9] (see also 

[ 7 ]). 


Theorem 2.1. Let (3) and (4) hold. Choose p = p{0,-) GR O RR(R;[0,1]) and 
y = y(0) = En[p]. Let y(-) be the solution of (2) with initial condition y. Define 
Pnit, ■) = Cn[y{t)]. Then pn converges almost everywhere and in Aoc([0^+oo) 
to the unique entropy solution p to the problem (5) with initial condition p(0, •) = p. 


3. Follow-the-leader model on net-works. In this section we extend the follow- 
the-leader model described in section 2 to a road network. We define a road network 
Q as a direct graph with W arcs (roads) and TVj nodes (junctions). We assume 
that for each junction J = 1,..., iVj , there exist disjoint subsets iNC(j), OUt(j) C 
{!,..., W}, representing, respectively, the incoming roads to J and the outgoing 
roads from j. Among junctions, we distinguish two particular subsets consisting of 
origins O, which are the junctions J such that iNC(j) = 0, and destinations V, which 
are the junctions J such that OUt(j) = 0. Finally, we denote by L^ the length of 
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the road R for any R = 1,..., TVr. At initial time, n vehicles are located anywhere 
in the network. 


3.1. A natural extension. A natural extension of the follow-the-leader model is 
derived as follows. We label the n vehicles by an index i = 1,... ,n and we denote 
by r(i, t) € {1,..., Nr} the road that the vehicle i is traveling at time t. We also 
denote by yi{t) the position of the vehicle i at time t, defined as the distance from 
the origin of the road 

There are three main differences with respect to the model on a single road, 
which require some modifications in the definition of 6i and w. 

A. The concept of “ahead” must be redefined because at junctions it is unclear 
which is the space in front of a car, unless a preference among the outgoing roads 
is assigned. We assume that each vehicle has a sequence of roads (i.e., a path) 
assigned at the initial time to travel along, so that it is always clear which is the 
space in front of the vehicle. We also assume that drivers can actually see the space 
ahead, so that they can always evaluate the distance between them and the car in 
front, even if the latter is in the next road. Given that, a vehicle is a leader if it is 
the first vehicle on its road and there is no other vehicles in the next roads along its 
path. It is important to note that now more than one leader can be present on the 
network at the same time, and vehicles can gain or loose the leadership. 

B. The vehicle i + 1 is no longer necessarily in front of vehicle i. This force us to 
introduce a new notation to refer to the vehicle in front. If a vehicle is not a leader, 
we denote by NEXT(t) the label of the car just in front of car i. If it is a leader, we 
set next(j) = 0. 

The distance between a (nonleader) vehicle and the vehicle in front is now com¬ 
puted as 


m ■■= 


2/NEXT(i)(i) if R{ht) = R(NEXT(f),t) 

(7^R(Lt) _ y .(t)) + Y,sL^ + J/nextW (t), if R(f, t) 7^ R(NEXT(f), t) 


where the summation with respect to S is done over all the empty roads between 
R(f,t) and R(NEXT(f),t) (if any) along the sequence of roads followed by vehicle i. 
It is important to note that Si(t) is in general discontinuous because at any time 
other cars can join or leave the i-th car’s path, thus reducing or increasing abruptly 
the distance 6i (or even changing the leader/follower status of the i-th car). 

C. The non-overlapping condition Si > in is no longer guaranteed, even if it holds 
at the initial time t = 0. To see this, let us consider the simple case of a merge, 
see Fig. 1. Assume that the outgoing road is empty and some cars are traveling 
along both the incoming roads (Fig. I-left). Assume also that the two leader cars 
are going to reach the junction nearly or exactly at the same moment (note that 
two cars are leader at this time). As soon as one of the two cars goes through the 
junction, the other one becomes a follower and its distance 5^, defined by (8), could 
be less than in or even 0 (Fig. I-right). However, two important observations can 
be made: 


Cl. The number of overlapping cars at a junction is bounded by the number of 
incoming roads of that junction. In other words, cars do not stack one over 
the other without limits, since the newly coming vehicles perceive the presence 
of the overlapping vehicles in front and behave normally, without getting too 
close to them. As a consequence, the total length of overlapping vehicles will 
tend to zero as n —?► -l-oo. 
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Figure 1. Overlapping issues across the junction. Leaders are 
denoted by filled square/circle. 

C2. The space interval in which overlap can occur is {L^ — which shrinks 

to 0 when n —> +oo. 

From now on we allow the above described “mild” overlap to happen and we handle 
it extending the function w defined in (3) by means of a new function w* defined 
as follows: 

w* : [0,+oo) [0,Umax], w*((5) = | (9) 

i.e. we assume that (smashed) vehicles with 5i < £„ stops completely until the 
vehicle in front leaves a space > t'n, then they re-start moving normally. If more 
than one vehicle is found to be exactly in the same place, the one with the largest 
label is taken as vehicle “in front” and the others “behind”. 

We are now ready to introduce the follow-the-leader model on networks. For any 
i = 1,..., n we write 

f i/iit) = if NEXT(i) ^ 0 

1 = ■Cmax, ifNEXT(i)=0 ^ 

with suitable admissible initial conditions yi(0) = yi. In addition, when yi{t) = 
^ the road R(*,t) must be updated according to the i-th vehicle’s path and 
yi{t) must be reset to 0 (similarly to a new initial condition). 

Remark 1. The extension of the follow-the-leader model described above defines 
implicitly a certain behavior of vehicles at junctions. Since vehicles cannot over¬ 
take each other, they basically adopt a FIFO behavior (see [22, Appendix B] for a 
general discussion in the context of traffic flow). Consider for example the case of a 
diverge (see Fig. 5(top-right)) and assume that one of the two outgoing roads is fully 
congested. If the first car on the incoming road wants to go to the congested road, 
it will stop. After that, all the following cars will stop too, even if they want to turn 
to the other outgoing road. The FIFO behavior will be kept in the macroscopic 
limit. 

Remark 2. Beside the overlap issue, the model exhibits other unrealistic behaviors. 
Let us consider the junction depicted in Fig. 5(bottom). If two cars arrive at the 
junction from distinct roads exactly at the same time and then they want to go to 
distinct roads, they will not be affected at all by each other, and, in particular, they 
will not slow down. This odd behavior can be easily avoided by considering enlarged 
(non point) junctions or allowing cars which approach a junction to see along all the 
incoming roads, but then it can be more difficult to prove the relationship between 
the microscopic model and the corresponding macroscopic model. 
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3.2. The model reformulated. In order to pass to the limit for n +oo, it 
is convenient to reformulate the model introduced in the previous section. Let us 
consider all the possible paths on Q joining all the origin nodes O with all the 
destination nodes T), see Fig. 2. Each path is considered as a single uninterrupted 
road, with no junctions (except for its origin and destination). Note that paths can 



Figure 2. A generic network Q. Two possible paths are high¬ 
lighted. Note that the two paths share an arc of the network. 

share some arcs of the network. Let us denote by P the total number of paths 
and let us divide the n vehicles in P populations, on the basis of the path they are 
following. 

Let us denote by Up, p = the number of vehicles following path p 

(we have and label univocally all vehicles by the multi-index (i,p), 

7 = 1,..., Up, p = I,..., P. Let us also denote by yi,p{t) the position of the vehicle 
i of population p at time t, defined as the distance from the origin of the path p 
(not from the origin of the current road). Since paths can overlap, a vehicle (fc, q) 
belonging to a population q ^ p with position yk,q(t) can be found along path p at 
time t. Let us denote by Cp(k,q,t) the distance between the vehicle {k,q) at time 
t and the origin of path p, along path p. In other words Cp acts as a change of 
coordinates between any path q ^ p intersecting path p and path p. The function 
Cp is trivially extended to vehicles already belonging to population p simply setting 

Cp{i,p,t) = y^^p{t). 

Most important, we redefine the concept of distance between two vehicles. First, 
we define NEXt(7,p; q) as the multi-index of the nearest vehicle belonging to popu¬ 
lation q in front of the vehicle {i,p) along the path p (setting it to 0 if there is no 
such a vehicle). Then, we define the p-distance between a (nonleader) vehicle and 
the vehicle in front as 


■ —?/NEXT(i,p;p) (^) 


This corresponds to consider only vehicles of population p, neglecting the presence 
of the others. We also keep considering the distance between two contiguous vehicles 
traveling along path p, regardless the population they belong to. To this end we 
define, for any nonleader vehicle of any population. 


:= Cp(NEXT(A:,9;*),t) - Cp{k,q,t), 


where NEXt(/c, g; •) is the vehicle in front of the vehicle {k, q), no matter the popu¬ 
lation it belongs to. 

In the new formulation the model (10) becomes. 



(11) 
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for any p = 1,..., P and i = 1,..., Up. Note that (11) is a system of P coupled 
systems of ODEs with discontinuous right-hand side. As already mentioned, the 
discontinuity comes from the fact that cars following paths other than p can abruptly 
join or leave path p, thus decreasing or increasing the distance However, in 

the time intervals in which no car crosses junctions the system falls in the standard 
theory of the ODEs. To deal with the discontinuity, we will consider the integral 
form of (11) and we assume that there exists a solution in the sense of Caratheodory. 

4. Micro-to-macro limit. To begin with, we note that condition (3) clarifies the 
relationship between the average density and the distance between vehicles which 
is valid in the limit n -|-oo (see Theorem 2.1). It states that 



( 12 ) 


In section 3.2 we have introduced p-distances which must be now related to the 
right average densities. To do that, let us use again (with an abuse of notation) the 
operators En and Cn introduced in section 2, modifying in a obvious manner the 
sets R (considering the total length £p of cars belonging to population p) and Yn 
to deal with the new framework. Define 


:= (yi.p, ■ ■ ■, 2/np,p), for any p = 1,..., P, 


and 



The new density functions are 





followers («,p) 
of pop. p 


p 





followers {k,q) 
on path p 
(any pop.) 


Both functions pO(t, •) and •) are defined on the whole pathp. The function 
corresponds to the density of vehicles belonging to population p, while the function 
uj^ represents the total density along path p, which depends on all vehicles belonging 
to population p plus all vehicles following a path q which shares some roads with 
path p, see Fig. 3. 

Dehning the limits (if they exist) 

fE(t,x) := lim p^{t,x), u!^{t,x) := lim ui'^(t,x), 

n^+oo n—)-+oo 

we get, in analogy to (12), 






(13) 


In order to introduce the macroscopic model we will also need the correspondence 
between the microscopic and the macroscopic velocity. To this end, we define the 
function v* such that 


V* : [0, -koo) —7- [0, U] 


'maxj 1 



(14) 
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Figure 3. Construction of functions (dashed blue) and ujP 
(solid red). Cars belonging to population p are marked by blue 
circles, all the others by red squares. Here £ = l,n = 9, np = 4. 


It is easy to see that the relation in (3) still holds, more precisely 

w*{S)=v* 

At this point it is crucial to note that a path, by definition, has no junctions 
(except for its origin and destination) and therefore is indistinguishable from a 
single road. As a consequence, we can follow in broad terms the results already 
proved in [7, 9, 29] for the micro-to-macro limit on a single road. 

Using (13), and reasoning by analogy with the ID case, we claim that, for any 
fixed p = 1,... ,P, the macroscopic equation associated to the model (11) is 

+ dx {pPv*{uj^)) = 0, t > 0, a; G R. (15) 

Note that in equation (15) both densities and appear. This is due to the fact 
that equation (11) describes the evolution of the population p only, but the velocity 
is evaluated considering the distance A^, which is indeed related to the total density 
on the path. Equation (15) accounts for all populations of cars and then it must be 
seen as a system of PDEs. 

To prove the correspondence, let us fictitiously extend each path to (—oo, +oo), 
setting the density to 0 after the leader and before the last follower. In this way we 
simplify the problem getting rid of the boundary conditions. 

Before entering the proof of the main result, we need to simplify the notations. 
Fix p G {1,..., P} and denote by i the index of a generic nonleader car belonging 
to the p-th population. Denote the index of the first car ahead of the i-th one, 
belonging to population p, simply by * + 1. Denote by k the index of a generic car 
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belonging to any population, including the p-th one. Denote the index of the car 
in front of the A:-th car by fc + 1. Finally, denote by ki the index of a generic car 
belonging to any population between car i (included) and car i + 1 (excluded), see 
Fig. 4(top); and by ik the index of the first car belonging to population p on the 
left of the car k {ik = k if the k-th car belongs to population p), see Fig. 4(bottom). 
Assume for simplicity that the leftmost car belongs to population p so that ik is 
well defined. 


i — 1 i ki i + 1 



ik k k + 1 *fc + l 



Figure 4. Pathp. Definition of ki (top) and definition of ik (bottom). 


Denote the initial density for the vehicles of population p by /2 G RCiBV (K; [0,1]) 
(dropping the superscript p for simplicity), and denote the corresponding initial 
positions of vehicles by {yi}i. Denote the positions of the Up vehicles of population 
p at time t by {yi{t)}i and the positions of vehicles of any population (p-th included) 
at time t by {yk{t)}k (note that their number is variable in time). Finally let us 
denote by T > 0 the final time for the simulation and consider a test function 
p € Cg°{{-oo,T] X R;R). 

We are now ready to prove the correspondence between the microscopic follow- 
the-leader model (11) and the macroscopic model (15), for any fixed p G {1,..., P}. 
Recalling [3, Formula (4.5)] the weak (integral) form of the equation (15) with the 
initial condition p,(0, •) = p, we compute 

pT 


fJ^nV*{uJn)(Px)dxdt-\- / fl{x)(f{0,x)dx = 
JR 

/ / V*{uJn)^x)dxdt x)ip{0, x)dx-\- 

Jo Jr ^ ^ Jr 


In := 

pT 


/ (fl{x) - fJ,n{0:X))(f{0,x)dx. 
Jr 

By definition, the densities p„ and uJn are constant in the interval [yk{t),yk+i{t)) 
for any t. Therefore we have 

pT D ryfc+i(i) / 




Pt + w 


Afc(t) 




dxdt + ^ 


rVi+i 


p(0, x)dx+ 


= Vk{t) 


(p(x) - p„(0, x))p(0, x)dx, 


where 5i = yi+i — yi- Denoting by p the total derivative of p w.r.t. time, we have, 
for x G [Pfc(t),yfc+i(t)), 

|pt(t,x) +yk{t)ipx{t,x) - p(t,yfc(t))| = 

Wt{t,x) + yk{t)px{t,x) - (pt{t,yk{t)) - yk{t)(pxit,ykit))\ < 

C(p{\ Umax) |yfc+l (t) yfc(0l 
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where the quantity Cip := Hv^Hc^ uniformly bounds from above the modulus of tp 
and all its derivatives up to second order. Equivalently, we can write 

Pt{t,x) + yk{t)px{t,x) = + 0{l){yk+iit) -yk{t)). 

Using the latter estimate in the expression for above and defining 

A ^ ^ f (p{0,x)dx, B := f - yLn{0,x))ip{0,x)dx, 


we get 

I^-A-B = ( Y. 


rVk+i{t) 


Kif) 


VkW 


^it,yk(t)))dxdt+ 


0 rVk+iit) 

0{l){yk+i{t) - yk{t))dxdt 

V Kit) U p) 


r YKt,ykit))p^,dt+o{iy^ rY^Kit)dt= 

Jo u (^) Jo ^ ^ik (0 


<1 


/ X] X] ‘Pit^ykiit))^j^dt + 0{l)in [ VAfc(t)dt = 

Jo Jq 


i fcig[i,j+l) 


X X {Kt^y^it)) + c>0)iy^+lit)-y^it)))^^dt+ 


i kie[i,i+l) 


o(i)C / X X Afci {t)dt 


^"X / (Kt^y^it)) + C>il){y^+lit)-y^it))) Y dt+ 


fcie[i.i+i) 


C>(l)^n X/ X Afc. (t)dt 

* fciG[*,i+l) 

iuY [ Kt,yrit))dt + 0{l)£r.Y [ Kt)dt + 0{l)£r.Y [ K^dt = 

I do ^ Jo j Jo 

^ Jo ^ Jo 

At this point we got rid of the populations g’s, p, and thus we are back to the 
case of a single road with a single population of vehicles. By the way, we have also 
resolved the discontinuity issues around the junctions. Therefore, from now on the 
proof follows the one in [7]. Recalling that <p(T, •) = 0, we have 

In = + 0{l)ln [ (^/rip (t) “ yi{t))dt AB = 

Jo 
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~^T I ‘p{0,yt)dx + O{l)en I ivn^it) - yi{t))dt+ 


p rVi+i 

I 





To proceed further, let us define A := — yi and note that it does not depend on 

Up OT n. We have 




Since y„(0, •) /2 (see [7, Prop. 2.1]) and in —t 0, all the terms in the latter quantity 

vanish as n ^ +oo. Moreover there exists y such that lim„^oox) = y(t,x) 
almost everywhere (see [9, Th. 3]). In conclusion, we have proved the following 

Theorem 4.1. Let (3) and (4) hold, and consider the extensions (9) and (14)- 
Fix p € {1,...,P}, and choose pP = fjP{0,-) £ R Cl iiy(R;[0,1]) and = 
En[p^]. Let y^(-) he the solution of (11) with initial condition yP (assuming the 
other solutions q ^ P are given). Define •) = C„[yP]. Then, pP{t,x) := 

lim„_>.+oo Pni'tjx) is a weak solution to (15) (assuming the other solutions qfi^p 
are given) with initial condition pP. 

Now we note that the system of PDEs (15) is nothing else that the multi-path 
model studied in [4, 5], 


dtp^ + dx{p^v(jjP)) = 0, p=l,...,P, t > 0, X € path p (16) 


where uiP(t,x) is the sum of all densities p‘^{t,x), q = 1,...,P, living at time t 
at the point x along path p and v is the velocity of the cars. It is a system of P 
conservation laws with discontinuous flux which provides an alternative method to 
deal with traffic flow on networks. It is characterized by the fact that junctions are 
embedded in the equations themselves, so that the dynamics at junctions have not to 
be resolved separately by ad hoc procedures, like, e.g., the maximization of the flux. 
It appears that the model maximizes the flux but under the additional constraints 
of equidistribution of the fluxes from the incoming roads (i.e. the incoming roads 
have assigned the same priority). 

Remark 3. We proved the correspondence in the limit between the microscopic 
and the macroscopic model population by population, i.e. Hxing p and assuming the 
other solutions to be given. This is different from considering the convergence of 
the fully coupled system of ODEs (11) to the system of PDEs (15). On the other 
hand, a complete study is much harder because the microscopic vector field is not 
regular and because theoretical results for the multi-path model are still missing. 


5. Numerical approximation. In this section we give some details about the 
discretization of the microscopic and macroscopic equations. 
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5.1. The follow-the-leader model. For numerical purposes, it is convenient to 
discretize the follow-the-leader model as it is described in section 3.1, rather than 
its equivalent formulation on paths. 

First we set Umax = 1 and 

y{p) = ^-P. w(5) = l-y. (17) 

Then, we introduce a time step At > 0 and we denote by y®, s = 0,1,2,... the 
approximate position y of the vehicle i at time t® = sAt. The discretization is 
obtained by the explicit Euler scheme 

y:+^ =y: + Atw*m^)), s = 0,1,2,... (18) 

for any car i = 1,..., n, and admissible initial condition y^ = yi. 

When a vehicle crosses the junction and goes beyond it, let say ?/i(t®) = l-|-e, 

it is assigned to the next road on the basis of its preference, and its position is 
updated as ?/i(t®) = e. We assume that no cars enter the road network after the 
initial time and that cars leave the network once they have reached a destination 
node. 

At the discrete level, a CFL-like condition is needed to ensure that vehicles do 
not bump into the one in front. Far from the junction, this condition is given by 

At w{6i) < 5i Vi. (19) 

Substituting (17) in (19) we get 

52 

At = At{n) < min --—. (20) 

i di- In 

Unfortunately, this condition does not prevent vehicle overlapping near the junction, 
as already discussed in section 3.1. In a discrete setting, the overlapping can happen 
also after the junction, more precisely it can happen either in the space interval 
{L^ — in,L^] of one incoming road R or in the space interval [0, At Umax) of one 
outgoing road. At Umax is the maximal distance one vehicle can travel in one time 
step. It can be seen that in the discrete setting the overlapping zone shrinks to 0 
for n —>■ -|-oo and At —>■ 0, and not only for n -boo. However, the overlapping 
is still “mild”, in the sense that the number of bumped vehicles is bounded by the 
number of incoming roads, see section 3.1. 


5.2. The multi-path model. For the numerical discretization of the multi-path 
system (15) we employ the Godunov-based scheme introduced in [4]. For the sake of 
completeness, we briefly recall the scheme. Along each path p, a space discretization 
is considered, defining a space step Ax > 0 and space nodes Xk ■= kAx, k = 
0,... ,NP, where is the number of space nodes on path p. We denote by 
and the approximations of the values yP{t’^,Xk) and uj'^it’^,Xk), respectively. 

The scheme reads, for any fc = 0,..., and s = 0,1, 2,..., 


s+l,p s.p 

Pk =Pk 


At 

Ax 


s.p 

Pk 

s,p 




0 S.p \ 

^‘^k+l) 


s.p 


UJ 


s,p 

k- 


1 


5(w 


s,p 

k-1 



( 21 ) 


with initial conditions and 


OJ 


s,p 

k 
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Note that is the sum of all densities living at Xk (along path p) at time t®. The 
function g is instead the classical Godunov flux 


9{P-iP+) 


min{/(p_),/(p+)}, if < p+ 

f{p-), if p- > p+ and p- < a 

f{<j), if p- > p+ and P- > cr > p+ 

f{p+), if p- > p+ and p +>(t 


where a := argmaXpgjQ f{p). Recalling (17), we have f{p) = p{l — p) and tr = i. 

We employ Dirichlet-zero boundary conditions /i^ = 0 at origin and destination 
nodes. 


Remark 4. The multi-path model, at least at the discrete level, does not allow the 
total density to be larger than 1 [5, Sect. 3.2]. This seems in contradiction with the 
possibility of overlapping in the microscopic model. Actually this is not, because 
the overlapping is “mild”, i.e. it concerns at most a finite number of cars, and then it 
is invisible at macroscopic level. However, a counterpart of the overlapping features 
at macroscopic level exists: the cells after the junctions act as a sort of buffer [5, 
Sect. 5.2], which is able to gather a flux of vehicles larger than maximal one, which 
is maxp{pu(p)}. The actual maximal flux is instead iVinc maxp{pu(p)}, where iVjnc 
is the number of incoming roads. 


6. Numerical tests. In this section we study the micro-to-macro limit by means 
of some numerical tests, finding a good agreement with theoretical findings of the 
previous sections. We consider the case of a simple network with 1 junction and (i) 
a merge, (ii) a diverge, and (hi) 2 incoming roads and 2 outgoing roads, see Fig. 5. 



Figure 5. Top-left: a merge (3 roads and 2 paths). Top-right: a 
diverge (3 roads and 2 paths). Bottom: a junction with 2 incoming 
and 2 outgoing roads (4 roads and 4 paths). 

In the cases (ii) and (iii), additional parameters are needed to describe the be¬ 
havior of the vehicles at the junction. They are usually referred to as distribution 
coefficients and specify the percentage of vehicles which wish to turn to the left and 
to the right. Since we are considering here simple networks with only 1 junction, 
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specifying the car’s turning conduct at junction is equivalent to define the car’s path 
on the network. In the case of more complex networks with more than 1 junction, 
we should instead assign the whole sequence of decisions each car makes at junc¬ 
tions. If we stick with per-junction coefficients, thus losing the global behavior of 
drivers, we expect convergence, in the limit, to the hybrid version of the multi-path 
model described in [5, Remark I]. 

Let us denote the distribution coefficients by Va^b, where a is the index of the 
incoming road and b is the index of the outgoing road. Clearly it is required that 

^Va^h = ^ Va. 

h 

In the follow-the-leader model, each vehicle traveling along road a is assigned to 
the road b with probability Va^b- In the multi-path model, these coefficients are 
used to define properly the partial densities ^^’s at the initial time given the initial 
condition for the total density. 

To improve readability and ease the comparison, we only show the total density, 
and we redefine it on single roads, rather than on paths. We also come back to the 
original notation p, denoting by the total density on each road R = I,..., ^Vr. 
We assume that all roads have the same length = 4 x 10^, and we divide them 
in = 100 space nodes. 

Comparison with the solution of the multi-path scheme is obtained by computing 
the average density of the microscopic vehicles, defined for each space node k and 
time step s as 

: y* e [a;fc,Xfc+i)}, (22) 

where 4hA denotes the cardinality of any set A. This quantity corresponds to the 
total length of cars found in each space cell divided by the length of the cell. 

6.1. Merge. In this section we consider a network with three roads and one junc¬ 
tion, with two incoming roads and one outgoing road. We denote by p^, p^ and p^ 
the density on the first incoming road, the second incoming road and the outgoing 
road, respectively (Fig. 5). Initial conditions are 

p^(0,x) = 0.5, p^(0, x) = 0.3, p^(0,x)=0 Vx. 

We run the simulation until the final time T = 3 x 10^. Fig. 6 shows the density 
computed by the multi-path model and the follow-the-leader model for two different 
choices of ^n and At. Since the outgoing road is not able to gather the fluxes coming 
from the incoming roads, two queues are formed and propagate back along the in¬ 
coming roads (with different speed). Incoming fluxes are equidistributed and queues 
have the same level of density, equal to p* such that p*v{p*) = ^av{a). Since we 
compare two numerically approximate densities, we expect a perfect match only for 
n ^ +O 0 and At, Ax 0. That said, numerical evidence confirms the convergence 
results, and also shows that the microscopic scheme is not diffusive as instead it is 
the macroscopic scheme. This is perfectly visible across the discontinuities. 

6.2. Diverge. In this section we consider a network with three roads and one 
junction, with one incoming road and two outgoing roads. We denote by p^, p^ 
and p^ the density on the incoming road, the first outgoing road and the second 
outgoing road, respectively (Fig. 5). Initial conditions are 

p^(0,x)=0.5, p^(0,x)=0, p'^(0,x)=0 Vx, 
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road 1 road 1 



Figure 6. Merge, result of the simulation at final time. Total 
macroscopic density redefined on roads (red line) and density of 
microscopic vehicles (blue circles). Left: £« = 3 (max number of 
cars per cell = 13.3), At = 3. Right: £„ = 1 (max number of cars 
per cell = 40), At = 0.2. 


road 1 


road 1 



Figure 7. Diverge, result of the simulation at final time. Total 
macroscopic density redefined on roads (red line) and density of 
microscopic vehicles T (blue circles). Left: in = 2 (max number of 
cars per cell = 20), At = 4. Right: in = 0.1 (max number of cars 
per cell = 400), At = 0.25. 


and distribution coefficients are 


= 0.8, lPi^4 = 0.2. 

We run the simulation until the final time T = 3 x 10^. Fig. 7 shows the density 
computed by the multi-path model and the follow-the-leader model for two different 
choices of in and At. As expected, the density splits among the two outgoing 
roads. The corresponding flux splits with ratio ^ and mass is conserved. Again, 
numerical evidence confirms the convergence results, however, the convergence is 
much slower than the previous case. This is probably due to the presence of the 
distribution coefficients V which introduce stochasticity in the system. Note that 
we are showing the outcome of a single run and not the average of many runs. 
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6.3. Junction with two incoming and two outgoing roads. In this section 
we consider a network with four roads and one junction, with two incoming roads 
and two outgoing roads. We denote by p^ and p^ the density on the hrst 

incoming road, the second incoming road, the first outgoing road, and the second 
outgoing road, respectiveiy (Fig. 5). Initiai conditions are 

p^(0,x) = 0.4, p^(0,x)=0.5, p^(0,x)=0, p^(0,x)=0 Vx, 

and distribution coefficients are 

iPl^3 = 0.7, iPi^4 = 0.3, iP2^3 = 0.6, iP2^4 = 0.4. 

We run the simuiation untii the finai time T = 4 x 10^. Fig. 8 shows the density 
computed by the muiti-path modei and the foiiow-the-ieader modei. Again, nu- 



road 2 



road 4 


1 

0.5 



Figure 8. 2-in-2 junction, result of the simulation at final time. 

Total macroscopic density redefined on roads (red line) and density 
of microscopic vehicles T (blue circles), in = 0.25 (max number of 
cars per cell = 160), At = 0.1. 

merical evidence confirms the convergence results, however, the convergence is even 
slower than the previous case. 

Conclusions and future work. In this paper we have investigated the many- 
particle limit for a natural extension of the follow-the-leader model on networks. 
Numerical tests have shown a slow convergence to the limit solution. This means 
that a large number of particles is required to get a precision comparable with that 
of the associate macroscopic model. In addition, the vector storing cars’ indices 
must be reordered after every time step in order to find the car in front of any car, 
so that a long computational time is needed. In conclusion, we discourage the use 
of the microsimulator but our results justify any multiscale approaches where the 
micro and the macro scale live together and exchange information. 

The convergence result proven in this paper is only one of the possible micro-to- 
macro limits which can be investigated in the framework of traffic flow models on 
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networks. First of all, one could define the proper follow-the-leader model which 
corresponds to, e.g., the LWR model with maximization of flux at junctions. This 
requires a nontrivial management of the junctions, in which an authority is able to 
decide who passes the junction and when. 

Second, one can consider analogous techniques for second-order models. In this 
case the problem is completely open because a second-order multi-path does not 
yet exist and then it can be difficult to find the right limit model. 

Finally, it could be interesting the investigation of the many-particle limit for 
meso-to-macro models. To this end, useful references could be the papers [11], 
which deals with kinetic models on networks, and [12] about fundamental diagrams 
for kinetic equations. 

Acknowledgments, The authors wish to thank F. S. Friuli, M. Falcone, and B. 
Piccoli for the useful discussions and suggestions. 
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